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ABSTRACT 

Electron degeneracy effects are dominant in ultra-dense plasmas (UDP), such 
as those found in white dwarfs. These effects can be treated systematically by 
obtaining an expansion of the screening length in inverse powers of h 2 . The theory 
exhibits Thomas-Fermi-like screening in an appropriate regime. In general, our 
theory leads to an 0(1) effect on the enhancement of fusion rates in white dwarfs. 
Further, it is shown analytically for these stellar conditions that Bose statistics of 
nuclei have a negligible effect on the screening length, in consonance with Monte 
Carlo simulations found in literature. 

Subject headings: screening, fusion, plasma, quantum corrections 
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1. Introduction 

The evolution of a white dwarf towards a black hole is marked by fusion reactions 
between carbon and oxygen nuclei. These reactions are enhanced by the screening of the 
Coulomb interaction between fusing nuclei by the surrounding plasma. Many papers have 
estimated this enhancement factor(Hamada and Salpeter 1961; Salpeter and Van Horn 1969; 
Janovici 1977; Slattery et al. 1982; Itoh et al. 1990; Ogata 1997; Ichimaru and Kitamura 1999; 
Pollock and Militzer 2004; Comparat et al. 2005; Gasques et al. 2005). It is difficult to gauge 
the accuracy of these calculations. Our goal here is to provide a systematic basis for the 
estimation of the fusion enhancement factor in white dwarfs. 

Nuclei in white dwarfs are submerged primarily in a sea of electrons. This is because the 
probability of finding a nucleus, say, in a two-component UDP (one specie of nuclei of charge 
Ze, and electrons) is down by a factor of 1/Z, compared to finding an electron. Furthermore, 
Coulomb repulsion between nuclei is stronger than the thermal energy, so that nuclei tend 
to remain separated over a relatively large distance. Consequently, we expect the screening 
to be dominated by a sea of electrons. These electrons are degenerate since the associated 
fermionic chemical potential is much larger than the thermal energy. 

Our earlier paper(Chitanvis 2007) was directed towards understanding the screening ef- 
fect of electronic quantum fluctuations on fusion reactions near the center of the sun(Salpeter 
1954; Gruzinov and Bahcall 1998). That paper showed quantum effects are negligible, via a 
systematic expansion of the screening length in powers of h 2 , putting to rest a controversy 
that has ebbed and flowed over the years(Bahcall et al. 2002). In this paper we consider a 
plasma where quantum effects dominate, which is quite opposite to the solar plasma. And 
so we contemplate the theory in our previous paper in inverse powers of h 2 , which will apply 
when the effects of electron degeneracy dominate. 

This paper provides an analytical underpinning to the numerical techniques used in 
the past to study fusion in UDP(Ogata 1997; Pollock and Militzer 2004). These papers 
show that the effect of indistinguishability of nuclei on the screening of fusion rates is small. 
Further, our technique provides an alternative, systematic method for the calculation of the 
enhancement factors of fusion rates in UDP. Our estimates of enhancement are much more 
conservative than those given by Ichimaru and Kitamura(Ichimaru and Kitamura 1999), and 
in general agreement with recent results of Gasques et al(Gasques et al. 2005). 
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2. Screening Formalism for degenerate electrons 

Let us begin with the classical Poisson-Boltzmann equation for a single species of ions 
and electrons: 

-V 2 = 4np 

P = P+ + P- 
p + = e n Z exp(—Ze(p/k B T) 

p_ = —e n Z exp(e0//csT) (1) 

where e is the magnitude of the electronic charge, ks is Boltzmann's constant, n is the 
average number density, Ze is the ionic charge, and T is the temperature of the system. We 
shall work in the linear regime, by retaining only terms first order in 0, leaving to the next 
section a discussion of nonlinear terms resulting from the Boltzmann distribution: 



= A 2 



A ° ~ *UirrJ(Z* + Z) (2) 

where A is the classical screening length. 

We shall now generalize this method to one where nuclei are treated classically, but elec- 
trons are treated quantum mechanically(Chitanvis 2007). The quantum-mechanical version 
of the linearized Poisson-Boltzmann equation for a single species of ions and electrons may 
be written in analogy with Eqn. 1: 



-V 2 = 4vrp 

P = P+ + P- 

( Zecj) 
P+ ~enZ[l- — 

P- =-em{r})\ 2 (3) 

where ip is the many-body quantum wave-function for electrons, and {r} refers collectively 
to the electrons in the system, and is the electrostatic potential. 
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We now invoke the following scaled variables, in order to ease subsequent calculations: 



4> = Ze<p/k B T 



A = ksT 



r 
A 



r' (4) 

Ak B T y ' 

where T' is defined differently from the usual plasma parameter. Notice also that the scalar 
potential has been scaled differently than in our previous paper (Chitanvis 2007). This is done 
to allow for a convenient analysis of higher order contributions in the next section. Note 
that the first of Eqns.4 shows that we are using /c B T as the energy scale. The electrostatic 
potential is then given by: 

V' 2 0= (0 + 47rr>| 2 -l) (5) 
This equation may be obtained from a Lagrangian density: 



C =-\\V4>\ 2 -v{^) 
v&fy =^ 2 + 4vr0r'H 2 -0 (6) 

The corresponding Hamiltonian density can be easily derived: 

H = ^\Vj>\ 2 + v(4>^) (7) 

We will now introduce second-quantized notation to deal with the statistics of electrons: 

v$, $) - v(4>, $ ± ) = - + 4vr0r'(^ + + (8) 

where ip± are Grassmann variables, and the subscripts refer to the spin of the electrons. 
The co-existence of Grassmann variables and scalars in Eqn.8 is not problematic, since we 
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shall use this discussion solely to define a partition function for the entire system. And soon 
thereafter we shall integrate over the electron degrees of freedom, so that only a functional 
involving the scalar potential survives. 

The total Hamiltonian H for the system, including the quantum-mechanical part for 
the electrons is: 



Ji = Ho + Hq 

H q = A (m)(|V^ + | 2 +|V^| 2 ) (9) 

The quantum correction has been encapsulated in the following dimensionless parame- 
ter: 

A ^1w) (10) 

where m is the mass of the electron, e its electronic charge, and ao is the Bohr radius. 

Since the temperatures we consider are T ~ O(10 8 K) ~ O(10keV), the mass density 
p ~ 10 9 g cm~ 3 , and the rest energy of the electron is 0.55MeV, it follows that the non- 
relativistic approximation employed in Eqn.9 is valid. 

The partition function may be written in scaled variables as: 

Z = J V0V 2 iP ± exp(- J d 3 x'(Ho + H Q )) (11) 
where it is understood that /c^T = 1 in the units we are using. 

Note that the total number of electrons associated with each ion of charge Ze is Z, and is 
obtained via < +t/>1'0_) >= Zn (where the angular brackets indicate an expectation 

value). We will indicate shortly how one may impose this constraint upon the system using 
a Lagrange multiplier. In condensed matter physics, this is done via the electronic chemical 
potential. In our problem, it will turn out to be more convenient to institute this constraint 
via a functional involving just the electrostatic potential. The quadratic nature of the energy 
functional in Eqn.ll allows us to perform the functional integration over the Grassmann 
variables associated with the electronic degrees of freedom(Ramond 1981), allowing us to 
obtain: 
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Z ~ J V$ exp(- J dV((l/2)|V0| 2 + (l/2)0 2 - $)) Det(JF) 

Det(.F) = exp(Trm(.T)) 

.F = — Ag(m) V 2 + AttT'4) (12) 

Having integrated over the electronic degrees of freedom, we are left with an effective 
energy density in terms of the electrostatic potential alone. We could have also done things 
the other way, integrating over the electrostatic potential in the partition function, leaving 
a quartic in the fermionic variables, as is conventionally done. Our procedure can be said to 
have bosonized our plasma, since we now only have the scalar potential to investigate. We 
shall show below how our method leads to useful insights into the statistics of the plasma. 

One way to impose charge conservation, which was discussed just below Eqn.ll, via a 
Lagrange multiplier by making the following addition (AH) to the energy density: 

AH = Aitvj) (13) 

Here v may be interpreted physically as a uniform charge density, which will be adjusted 
to ensure either overall charge neutrality. We now need to evaluate the determinant of 
the operator obtained in the process of performing the quadratic functional integral over 
fermionic variables. This is conveniently performed in Fourier space. 

/H 3 k 
— ln(4vrr'0(A;) + A Q (m) k 2 ) (14) 

where <f>(k) is the Fourier transform of <p. 

For temperatures in the range of 2 x 10 8 i^, and using a number density for our plasma 
~ 2 x 10 32 cm -3 (which corresponds to a mass density of 3 x 10 9 g — cm" 3 ), we find with 
Z ~ 10, it turns out that A ~ 1 Fermi(10- 12 cm), V ~ 60, and A Q (m) ~ 0.3 x 10 5 . This 
is an indication that quantum corrections are dominant in this system. Hence AuT'/Aq ~ 
2 x 10~ 2 << 1 is an excellent choice for a perturbation parameter: 

ln(4vrr'0(A;) + A Q k 2 ) « ln(A Q k 2 ) + - ( + 0((A.T> / A Q f) (15) 

Q \ Q J 

Terms devoid of field variables in Eqn. 15 will be ignored. For the particular case of 
a dense UDP considered here, the cubic term ignored in the above expansion is of O(l0~ 2 ) 
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compared to the bilinear terms we retained. Thus, higher order terms can be safely ignored. 
We must also consider for consistency higher order terms in the expansion of the Boltzmann 
distribution corresponding to the nuclear charges. We shall do so in the next section. 

Note the appearance of powers of k~ 2 in Eqn.15. This could be problematic when we go 
to a representation in real space, as Fourier Transforms of k~ 2p , with p — 2, 3, ... do not exist. 
For this reason, we make the following ansatz for <p(k) which resolves this issue entirely: 



4>(k) = k 2 X (k) 

0(f) = -V\(f) (16) 



-2 



Then the corresponding Lagrangian density in real space may be written to order A Q 
as follows: 



^effective = ~ ~ | VV 2 x(r) | 2 

-(1 + 47n,)V 2 x(f) + i(V\(f)) 2 - V effective ( X (^) 
VeffectUx^)) « bxlf) + b 2 xif) 2 (17) 

where b = AttT/ Aq. v will be chosen to zero out the coefficient of the term linear in x, in 
order to maintain overall charge neutrality. The equation of motion obtained by extremizing 
the above Lagrangian with respect to variations in the field x is linear, albeit of order six: 



(A 3 - A + b 2 ) X {r) =-b 

A = V 2 (18) 

The particular solution of the inhomogeneous solution of Eqn.18 turns out to be, by 
inspection, simply — b^ 1 . This constant is physically meaningless, as it does not contribute 
to the scalar potential and hence to the charge density. To this particular solution must be 
added the solution of the homogeneous part of the differential equation. Upon factorizing 
the trinomial above, the homogeneous part of the differential equation may be cast as: 



(A - s 3 )(A - s 2 )(A - s 1 ) X (f f ) = (19) 

where Si,s 2 ,S3 are the roots of the trinomial in Eqn.18. These can be found can be done 
most easily using Mathematica. Upon careful examination of these roots, only one yields the 
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correct limit for the screening length given by l/y/s 1 , l/y/s 2 , 1/ y/s z m ^ ne limit that K — > oo. 
The screening length must go to infinity in this limit due to quantum fluctuations.- One 
must remember that quantum fluctuations increase the screening length(Chitanvis 2007). 
We will denote this root by s\. Its full expression is uninformative. However, in limiting 
cases it may be written down simply. 

For H — > oo, 

si « b « 1 (20) 

For this limiting case, the dimensionless screening length is then l/y/b >> 1. 
Hence we shall choose the physically interesting solution to satisfy: 



(V 2 - Sl ) X (r) = (21) 

and this will be sufficient to guarantee that the sixth order differential equation is auto- 
matically satisfied. One may verify the screening length argument by noting that in one 
dimension exp(— x/y/b) solves the differential equation in the limit that h — > oo. We have 
used here a finite-temperature formalism involving not only an electron gas, but also an 
admixture of positive nuclei to study screening in the limit that electronic quantum fluctu- 
ations dominate. We therefore expect to get for h — > oo a screening length which is similar, 
but not identical to the Thomas-Fermi length, which is obtained from a zero-temperature 
formalism. As given in Fetter and Walecka(Fetter and Walecka 1971), the Thomas Fermi 
screening length for an electron gas is given by: 



I Cipro 

^Thomax- i '< rim ' s "\ / 2 

h 2 e 2 

"° (22) 

where do is the Bohr radius, and r is the mean free distance between charges. 
In our case, the screening length 1/y/b (h — > oo) turns out to be: 



^ ..,,>-o(T) 

/v screening ~ W ^ u^. ;t 

£ (T) = ^ (23) 
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where the temperature-dependent scale £o(T) plays the same role as does r in Eqn.22. 
Notice that the screening length displays a rather weak ~ T 1//4 temperature dependence, 
in qualitative agreement with the conclusions of Ichimaru(Ichimaru and Kitamura 1999) for 
electronic screening in UDP. In general, the screening length is given by l/y/s 1 and will be 
discussed in much detail in of section 4. 



3. Non-linear screening effects of nuclear charges 

In the previous section, we discussed in detail how non-linear terms arising from electron 
degeneracy are safely of much smaller magnitude than those retained. The issue appears 
to get turned around when one considers non-linear terms arising from higher order terms 
in the expansion of the Boltzmann distribution corresponding to nuclear charges. This is 
because the coefficients of the nonlinear terms are larger than those found for degenerate 
electrons in the previous section. The characteristic Coulomb energy Z 2 e 2 /r >> kT, for 
Z = 6. Hence we are in the regime of strong screening(Salpeter 1954). What this means 
is that the positive charges remain separated by a large distance, with electrons performing 
most of the screening of nuclear charges. We will see shortly how this is borne out somewhat 
more rigorously. 

The issue may be compactly discussed by noting that retention of the quadratic term 
in Eqn. 3 leads to the following cubic modification of Eqn.8: 

v v - ^ 3 (24) 

The corresponding partition function for our classical system resembles one for a Eu- 
clidean scalar quantum field theory (Ramond 1981). Effects of the non-linear terms can be 
evaluated perturbatively, using Feynman diagrams. Feynman diagrams can be used to esti- 
mate the leading order contribution from this cubic term to the self-energy of the system (or, 
equivalently, the dielectric constant, or the screening length). It turns out that the actual 
value of the contribution is numerically quite small. This suggests that the linear screening 
approximation we retained in the previous section is an acceptable approximation. 

In order to perform this calculation, we formally ignored terms of 0(b) ~ 10 -2 << (1/6) 
in our diagrammatics. We then used Mathematica to formulate symbolically the lowest order 
contributions from the cubic potential. The term that survives is a polarization-like diagram 
which comes from terms of (9((l/6) 2 ). The calculation is done as usual in momentum space. 
At this point the momentum variable is simply set to zero, so that we get an expression for 
the inverse square of the screening length: 
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^(k) p olarization-like(k) — — J j^^ Go(k)Go(p — fc) 

Go(p) = ^ (25) 

where S(/c) is the usual self-energy. It can be easily related to the dielectric constant of the 
UDP. We shall restrict attention to k — 0, when the self energy it reduces to the inverse 
square of the screening length, and is sufficient to allow us to gauge its magnitude relative 
to the degenerate contribution in the previous section. 



Whence 

^-polarization-like = = 0) polarization-like 



967T 

~ 3.3 x 10~ 3 (white dwarfs) 
^ P oiarization-uke ~ 1.7 x 10~ n cm (white dwarfs) (26) 

We see that the screening length from the cubic term in the energy functional (due to 
nuclear charges alone) is larger by an order of magnitude than that obtained in the previous 
section by an order of magnitude (A screening ~ 10~ 12 cm for a white dwarf, from Eqn.19). 
This difference in the screening due to nuclear and electronic charges means that electrons 
can be in closer proximity to a nucleus than another another nucleus. 

The screening lengths appears as its inverse square in the corresponding differential 
equation for the scalar potential. Hence one can combine the contributions from the previous 
section, with the current one: 

^-polarization— like ^screening) '^(0 ^ (27) 

Thus the lowest order contribution of the cubic term to the screening length in Eqn.27 
is about two orders of magnitude smaller, compared to the screening contribution obtained 
in the previous section. We shall therefore ignore this non-linear correction. 

The conclusions of this section may have to be modified if higher order terms turn out to 
be larger. For that purpose, we retain the quartic term in the energy functional arising from 
the Boltzmann distribution for the nuclear charges, so that Eqn.24 is modified as follows: 

« - v - + (28) 
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The lowest order contribution from the quartic term is the setting-sun diagram (Ramond 
1981). There are no cross-terms at this order between the cubic and quartic terms. 

1 C f cftk (ftk - - -> -> 

S setting- S un(p) = g J Jjj^jj^G (k 1 )Go(k 2 )Go(\p-k 1 -k2\) (29) 

Using t -1 = f d\exp(—\t), converting the momentum integrals to center of mass and 
relative co-ordinates, and using dimensional regularization, we can compute the self-energy 
for p = in the following form: 



^ setting — sun ( ) 



1\ L I 



6 8tt 4 



x L ~ e 



i - t HO-Hi + y) dy (30) 



where £ — > + , a small-distance cut-off, and e — > + have been inserted to guarantee conver- 
gence. Using Mathematica, it can be shown that: 



L = K)-T( £ ) 

7 = hjfl _ 1 / 3- 4 cosec( Te ) + jFi(1| 1; 2 _ £ , p _ (:!;| ) 



4e 16 

where 2-^1 (a, c; z) is the hypergeometric function. 

Employing Laurent expansions within Mathematica in powers of e, and upon using 
counter-terms that account for divergences, one obtains the following correction to the square 
of the screening length, correct to lowest order in the quartic term: 



^setting— stin(O) 

^ (-7 - MO + Is (" 2 * 2 + (5)) ~ Ti (( 21og2(3) " -»») 

V 7 ? = I (32) 

where a « A is a microscopic length cut-off required to render the integrals finite, 7 ~ 0.577 
is Euler's constant, and ^1 « 0.572 is the value of the derivative of the Hypergeometric 
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function which appears in Eqn.31. A logarithmic dependence of our answer on this cut-off 
implies a relative insensitivity to this parameter. It is clear that the theory used in this 
section is certainly not valid at the nuclear level, and so we will use £ = a/A ~ 10 2 . A 
cursory examination shows that this particular diagram yields a small contribution (9(10~ 4 ) 
contribution to the self-energy. Then following the argument above, for the cubic term, this 
screening contribution can be ignored as well. Thus there appears to be a trend for the 
classical, nonlinear terms to be small. 

Of course, further issues may arise in this perturbative argument as even higher order 
terms arising from the Boltzmann distribution are contemplated. We will leave these ques- 
tions for future investigation. It may be possible to extend to white dwarfs the methods 
utilized by Brown et al(Brown et al. 2006) for obtaining the screening length for a dilute, 
highly charged plasma. 

4. Quantum effects of nuclei on screening 

Over the years, researchers have delved into the importance of applying a quantum- 
statistical treatment to the nuclei surrounding the ones undergoing fusion. This is a reason- 
able point to investigate, given the extremely high densities available in white dwarfs. Path 
integral Monte Carlo (PIMC) techniques have led to the discovery that the bosonic nature 
of nuclei(Itoh et al. 1990; Ogata 1997; Pollock and Militzer 2004) makes a calculable, small 
contribution to the screening length, or equivalently, to the dielectric constant of the UDP 
under consideration. Here we provide an analytical underpinning to that observation. The 
second point that needs to be reinforced is that quantum effects of nuclei in the UDP are 
small in general, compared to the electronic contribution. The argument is basically that 
of Born and Oppenheimer, who showed that the nuclear mass is so large compared to the 
mass of an electron that an adiabatic approximation can be applied. That argument has to 
be extended to finite temperatures. 

We begin by considering a UDP consisting of spin-zero nuclei, in addition to a sea of 
degenerate electrons. Thus the Poisson-Boltzmann equation (Eqn.8) will be replaced by: 

v(j>, -> v B (4>, *, t/>±) = 47rZ0r(^ f *) + 47r0r<>^ + + $L-0_) (33) 

where ^ are the second-quantized creation and annihilation operators corresponding to 
bosonic nuclei, and we have continued to assume that the photons are numerous and hot 
that they can be treated classically. 

Equation 9 must be modified to account for the free Hamiltonian of the nuclei: 
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Aq(M) |W 



(34) 



where Aq(M) is defined via Eqn.10, with the mass of the electron being substituted by the 
mass of the nucleus under consideration. It has been assumed that the interaction between 
nuclei are dominated by the Coulomb potential. This is a reasonable assumption to make in 
view of the fact that in this dense plasma, the nuclei are separated on the average by about 
lOFermis, as shown in the previous section. 

The corresponding functional integrals involving the nuclear field variables can be done, 
just as the fermionic degrees of freedom were accounted for. The net result for the partition 
function is: 



Using parameters representative of a white dwarf, and given that M ~ 2 x 10 4 m, Z ~ 10, it 
turns out that A Q (M) ~ 10~ 4 A Q (m), so that A Q (M)/4nZT ~ 1(T 3 . This clearly shows that 
quantum effects, including Bose statistics of nuclei in the UDP are negligible in circumstances 
representative of a white dwarf. Thus, we may continue to use a classical treatment for 
nuclei, as was done in the previous section. One can regain from Eqn.36 the Boltzmann 
approximation using the method outlined by Chitanvis(Chitanvis 2007). This will yield 
systematic, extremely small quantum corrections to the Boltzmann approximation in powers 
of h 2 . Considering the nuclei as fermions does not change this conclusion. 

Quantum effects of nuclei are much smaller than the quantum effects due to electrons, 
primarily due to the large mass difference. As such they can be ignored. Our conclusions 
concerning the effects of indistinguishability are in general agreement with Itoh(Itoh et al. 
1990), Ogata(Ogata 1997) and Pollock and Militzer (Pollock and Militzer 2004). These au- 
thors did not consider the effects of electron degeneracy on the same footing as the nuclear 
quantum effects. 




(35) 



where 




(36) 
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5. Results 

We shall compare our results with those found in published literature. For succinctness, 
we shall quote the enhancement factors of pycnonuclear reaction rates given in the review 
article of Ichimaru and Kitamura(Ichimaru and Kitamura 1999) and results by Gasques et 
al(Gasques et al. 2005). 

Our approach provides an integrated, first-principles theory of screening effects due to 
electrons and associated nuclei. Furthermore, we have derived a systematic way of obtaining 
screening effects in a UDP. As such, it is possible for us to estimate the accuracy with which 
we can calculate our enhancement factors. In order to do that for white dwarfs, we must 
first generalize our formalism to a binary ionic mixture (BIM), e.g. a mixture of 12 C, 16 O 
nuclei and associated electrons, as representative of the species in a white dwarf. This is 
easily accomplished via the following substitutions in all our formulae: 



z 2 _, z 2 = niZ * + n2 ^ 

n 

n — > n = Hi + n 2 (37) 

where rii, n 2 are the average number densities of each species in the UDP. These substitutions 
arise naturally through a re-derivation of our theory for a BIM. Generalizations to more than 
two components is straightforward. 

Ichimaru uses improvements over a standard procedure to obtain net enhancement fac- 
tors. First, the screening length/dielectric constant of an electron gas in a jellium of positive 
ions is obtained. Then various sophisticated methods are used separately to obtain screening 
effects due to nuclei surrounding the moieties undergoing fusion. Physically reasonable mix- 
ture rules are utilized to obtain the overall enhancement of nuclear rates in a UDP, caused 
by the screening of the nuclear Coulomb repulsion by intervening charges. It is not possible 
to gauge the accuracy of such calculations. In general we find our estimates for fusion rate 
enhancement in white dwarfs are much more conservative than those of Ichimaru and Ki- 
tamura (Ichimaru and Kitamura 1999). In an effort to understand this difference, we found 
that the overall plasma screening length obtained by us (primarily a degenerate electron con- 
tribution), and the electronic screening length of Ichimaru and Kitamura agree reasonably 
well. It follows that the difference in the enhancement rates must come from their treatment 
of screening due to the nuclei. And we showed in section III that our classical treatment 
of the heavy nuclei gives a very accurate result for the screening correction. Furthermore, 
the values of the classical plasma parameter T dassica i = Z 2 e 2 /r kT « 170 (r is the mean- 
free distance between particles) for the cases listed in Table 1. As such there is no concern 
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regarding the UDP being close to a crystallized state. 

Gasques et al utilize a re-parameterized version of the enhancement factor obtained 
by Slattery et al(Slattery et al. 1982). Our results agree more reasonably with Gasques et 
al(Gasques et al. 2005) for the following reason. The enhancement factor of Gasques et 
al is defined in terms of the plasma parameter of a one-compenent plasma. In our case, 
we have a binary mixture. However, since we are in a strong-screening regime(Salpeter 
1954), we suggest that our BIM is a one-component plasma composed only of electrons (for 
purposes of screening). In any event, it is instructive to compare the enhancement using a 
classical electron gas to our results. We find that the classical result gives an enhancement 
factor which is roughly twice that obtained for our BIM. We attribute our lower values of 
enhancement to quantum effects. 

An interesting observation to be made is that if we use the enhancement formula of 
Gasques et al(Gasques et al. 2005) for a one-component classical plasma composed of positive 
charges, each having a charge Ze, with Z = 6, there is remarkable agreement with the results 
of Ichimaru and Kitamura(Ichimaru and Kitamura 1999). 

The numerical values obtained for our screening length and comparisons to Ichimaru's 
enhancement factor(Ichimaru and Kitamura 1999) and that of a classical electron gas(Gasques 
et al. 2005) have been encapsulated for white dwarfs near ignition in Table 1. 
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Quantity 


Current theory 


(Ichimaru and Kitamura 1999) 


classical electron gas 


Model 


WD1 






Density 
T (Kelvin) 
Composition 


<- 


3.0 x 10 9 g - cm~ 3 
1.8 x 10 8 K 

12 C, e~ 


-> 


^enhance 


0.57 


12.09 


0.26 


Model 


WD2 






Density 
T (Kelvin) 
Composition 


<- 


9.0 x 10 9 g - cm~ 3 

1.1 x 10 8 K 

12 C, e~ 


-> 


■p 

1 enhance 


1.39 


23.10 


0.73 


Model 


WD3 






Density 
T (Kelvin) 
Composition 




9.0 x 10 9 g - cm" 3 
3.4 x lO 7 ^ 
12 C, e~ 




"P 

1 enhance 


6.02 


20.76 


2.5 


Model 


WD4 






Density 
T (Kelvin) 
Composition 




9.0 x 10 9 # - cm~ 3 

1.1 x lO 8 ^ 

12 C (75%), 16 (25%)e" 




^enhance 


1.42 


23.12 


0.71 



Table 1: Comparison of quantum corrected rate factors for white dwarfs near ignition, be- 
tween our calculation and two previous calculations(Ichimaru and Kitamura 1999; Gasques 
et al. 2005). The enhancement of the fusion rate is calculated as exp(r en / iance ), T en h ance = 
e 2 / (A screening kBT) . The quantum-influenced screening length A screening is defined in dimen- 
sionless terms as via Eqn.19. The different scenarios for white dwarfs, viz., models 

WD1 — WD4 are described in Ichimaru's paper(Ichimaru and Kitamura 1999). 



